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Abstract 

We review some occurrences of Painleve II transcendents in the study of two-dimensional 
Yang-Mills theory, fluctuation formulas for growth models, and as distribution func- 
tions within random matrix theory. We first discuss settings in which the parameter 
t£| ■ a in the Painleve equation is zero, and the boundary condition is that of the Hasting- 

MacLeod solution. As well as expressions involving the Painleve transcendent itself, 
one encounters the sigma form of the Painleve II equation, and Lax pair equations 
in which the Painleve transcendent occurs as coefficients. We then consider settings 
which give rise to general a Painleve II transcendents. In a particular random matrix 
setting, new results for the corresponding boundary conditions in the cases a = ±1/2, 
1 and 2 are presented. 

1 Introduction 

A Painleve II transcendent refers to a solution of the Painleve II nonlinear differential 
equation 

q" a = tq a + 2q 3 a - a. (1.1) 

For special values of the parameter a, in particular a = 0, and special boundary conditions 
determining a particular solution of (II. ip . Painleve II transcendents are prevalent in random 
matrix theory. They come about in the so-called double scaling limit, taken in the study of 
matrix models (see e.g. [Sj), as the scaling limit of certain probabilities specified as averages 
over unitary, orthogonal and unitary symplectic ensembles, as well as the scaling limit of 
the gap probabilities for ensembles of Hermitian matrices at the soft edge, when the leading 
order eigenvalue density goes to zero like a square root (these latter two occurrences will 
be reviewed below). As also to be reviewed below, it is furthermore the case that certain 
Lax pairs or Painleve II transcendents arise in the context of the study of the scaling limit 
of some generalizations of the above settings which involve an additional parameter. 

Our primary interest in this paper is to highlight the occurrence in random matrix 
theory of (11.11) for general a. However, we take the opportunity to also review more 
generally specific settings in random matrix theory and its applications leading to Painleve 
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II transcendents. In this regard, low- dimensional field theories, which give rise to random 
matrix (or random matrix-like) averages, will be considered first. These very same averages 
appear in the study of some probabilities associated with certain combinatorial models and 
growth processes, so these will be considered next in Section 2.3. In Section 2.4 we turn 
to the problem of soft edge gap probabilities in Hermitian random matrix ensembles with 
unitary, orthogonal and symplectic symmetry. The corresponding formulas in terms of 
a Painleve II transcendent for the probability of no eigenvalue in (s, oo) are precisely the 
same as those appearing as limiting probabilities of the combinatorial models from Sections 
2.2 and 2.3. An explanation for this involving first a transformation identity between 
probabilities, then an analysis of the hard-to-soft edge transition in random matrix theory 
as given in [13] is revised. In Section 2.5 we describe settings, within combinatorial and 
growth models, and within random matrix theory at the soft edge, which give rise to Lax 
pairs for a Painleve II transcendent. 

The Painleve II transcendents in the above problems either literally satisfy the case 
a = of (jl.ip . or can be mapped to that case. In Section 3 we turn our attention to a 
setting in random matrix theory at the soft edge in which a is a continuous variable. After 
presenting some review material we proceed to give new results in the cases a = ±1/2, 1 
and 2. 

2 A brief survey of Painleve II with a = in random 
matrix theory and applications 

2.1 Two-dimensional Yang-Mills gauge theory 

In two- dimensions QCD without quarks is a pure Yang-Mills gauge theory. The setting 
involves a two-dimensional orientable manifold M. of volume form y/g, containing at each 
point x a gauge field (Ai(x), ^(x)). The latter are taken to be N x N matrices from the 
Lie algebra of an appropriate gauge group G, and from these field strengths are defined by 



A fundamental property of the field strengths is that under the local gauge transformation 



they transform according to F^ u i— > S 1 (x) F^ v S \x) . This in turn demonstrates the gauge 
invariance of the Yang-Mills action 



F, 



flU 



dfiA v — d v A^ + z[v4^, A u ] 



(/i,z/=l,2). 



(2.1) 



A^ ^ S ,_1 (f)A /1 5'(f) - iS'^x^Six), S(x) e G, 




where 1/A 2 is a coupling constant, and also of the partition function 




(2.2) 
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The notation G in (12. 2p denotes the particular gauge group under consideration, while 
[XM.^] denotes a suitable measure on the fields. 

One approach to the computation of (12. 2 p is to discretize the manifold and also the 
action [19]. In two-dimensions the manifold can be triangulated into basic units called 
plaquettes. On each edge L of the triangulation, and with gauge group G, one associates 
an element Gl G G and defines the corresponding lattice regularized partition function 
according to 

z jttto = J Y[{dG L ) J] Zp(G P ), (2.3) 

L plaquettes 

where (dG) denotes the Haar measure on G. Here Gp = rLepiaqucttcs is the Wilson 
loop product [68]. Taking Zp as the Wilson action, 

Z P (G P ) = exp (bNTr(G P + G P )j , (2.4) 

which in the case that G is a classical group is a valid choice since it reduces to e~ lYM in 
the continuum limit, one sees that (12. 3p factorizes as 

Z jattice = n e WTr(G L+ Gl)\ \ M (2.5) 

\ \ / GlGG/ 

where M is the number of plaquettes. Hence we have the result that (^k^ tlce )V M i s given 
by an explicit matrix integral. We will revise below how a certain double scaling limit - 
taking N — > oo but also tuning b as a function of N — leads to a Painleve II transcendent. 
In the context of the present problem the double scaling limit is motivated by the desire 
to relate two-dimensional lattice QCD to string theory; see e.g. [55] . 

The choice of action (12.41) is not unique if the only requirement is that it reduces to 
6 -Iym j n the continuum limit. In addition to requiring that Zp[Gp) reduces to the Yang- 
Mills action, following [415] one can also require that fusing two adjacent plaquettes P\ and 
P2, and integrating out the group element associated with the common edge give back the 
action Zp 1+ p 2J 

J \dGs) Zp 1 {GiG2G^)Zp 2 {GiG^ ) G\) = Zp 1+ p 2 {GiG2GiG^). (2.6) 

This has the unique solution Zp{G) = (G|e~( Ap / 2Ar)A |II), where Ap is the area of the 
plaquette P and A is the Laplacian for the classical group G. Moreover, expanding about 
U 1 — y I -\- iA^dx^ 1 shows that in the continuum limit Zp reduces to e~ lYM . 

The property (12. 6p is used to argue that the corresponding lattice representation is in 
fact exact for any triangulation. Using the smallest number of triangles consistent with 
the topology of Ai, Z M can thus be computed exactly. Consider in particular the case 
that M. is a sphere of area A. One first notes that this quantity together with the rank N 
of the classical group can be taken as the parameters in (I2.2p . allowing the coupling A 2 to 
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be scaled out of the problem by setting A 2 iV = 1. Writing then Z M (\; G) = Z(A; G), and 
with 

A(xi, . . . ,x N ) = (x k -xj), (2.7) 

l<j<k<N 

the Vandermonde product, the exact partition functions obtained this way for the group 
G equal to one of the classical groups are [611 [T9] 

oo 

Z(A;U(N)) oce-^ 2 - 1 ^ 24 ^ A 2 (m, . . . , n N )e^ A ' 2N ^=^ 

ni,...,njv=— oo 
oo N 

Z(A; Sp(2JV)) oc e ^+V2)(iv + i)/i2 ^ ^ ^ )( "Q n 2 )e -{A/w) E f =1 „j 

ni,...,n JV =-oo j=l 

oo 

Z(A; + (2iV)) oc e ^-i/2)(AT-i)/i2 ^ A 2 (n 2 ; . ^^g-W^Ef^n^ (2 _ g) 

n± ,. . .,njy=— oo 

where the proportionality constants depend only on JV. Like the matrix integrals in (12.51) . 
these expressions can be analyzed in the double scaling limit and give rise to a Painleve 
II transcendent. In fact the functional forms obtained are identical to those for the corre- 
sponding case of (I2.5p . as will now be reviewed. 

2.2 The double scaling limit 

Consider for defmiteness the group integral in (12.51) in the case G = U(N). In the limit 
JV — > oo this exhibits the leading order form [H] 

ton _Llog/e 6 ^(^ + )\ = J 3 1 oL L 1 2 (2.9) 

n^ooN 2 6 \ /t/ e u(iv) I 26 log26, 6>± v ; 

k 4 2 2 

and so is discontinuous in the third derivative at 6 = 1/2. The first rigorous proof of 
this result was given by Johansson [44]. As mentioned, a string theory viewpoint of two- 
dimensional QCD motivated tuning the coupling 6 to scale with JV as it approaches 1/2. 
Specifically, setting 6=1/2- 2- 4 / 3 JV" 2 / 3 t it was shown by Periwal and Shevitz [57] using 
orthogonal polynomial methods that then 

^ lim ( log (e^^) - 6 2 JV 2 ) = -g 2 (t), (2.10) 

where qo(t) satisfies the Painleve II equation (II. ip in the case a = 0. Of course one still has 
to determine the boundary condition. This was done by Crnkovic et al. [20], who showed 
that it is required q (t) — > (— \t) l l 2 as t — > — oo and q (t) — > as t — > oo. Moreover, it was 
pointed out in [20] that the existence and uniqueness of such a solution of Painleve II has 
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been known since the work of Hastings and Macleod [12] . In the latter work the condition 
q (t) — > as t — )■ oo is further refined to 

q (t) ~ Ai(t), (2.11) 

t— >oo 

where Ai(£) denotes the Airy function. 

Consider next Z(A;\](N)) as given in (12. SJ) . It was shown by Douglas and Kazakov 
[21] that 

fe^w; u w>={£$; 115 < 2 - i2 > 

The function F_ can be expressed explicitly in terms of elementary functions, while F + has 
an explicit form in terms of elliptic functions. For present purposes, the important feature 
is that when expanded about A = tt 2 the two functions agree up to order (A — vr 2 ) 2 , and 
thus like (12. 9 j) there is a third order phase transition. Gross and Matytsin [3D] analyzed 
the corresponding double scaling limit, obtaining the result 

\im ^\oge- A2N2 ^Z(A-MN)) . . = -«?(*), (2-13) 



J 4=7r 2 -f7r 2 /(27V) 2 / 3 



where go(^) is the very same Painleve II transcendent as in (12.101) . 

The double scaling limit of the group integral in (12. 5 p in the cases of the orthogonal, or 
unitary symplectic classical groups was first undertaken by Myers and Periwal [52]. They 
obtained the results 



and 

lim ( log (e bNT *°) - 2b 2 N 2 ) = lim ( log (e bN ^ s ) - 2b 2 N 2 

N^oo dt 2 \ \ /OeO+(27V+l) / N^oo dt 2 \ \ I 5eSp(2Af) 

1 



where b is scaled as specified above ( 12 .101) but with N replaced by 2N. We observe that 
adding ( 12.14ft and ( I2.15jl gives (12.10jl . This can be understood as a result of the more 
general identity [6] 

2n+l n n+1 

( " i < ei9j )) = (n°( e "o) ~ (ii< eiej )) — ( 2 - 16 ) 

\ AA / U(2n+1) \- LJ - /Sp(2n)\- LJ - /0+(2n+2) 

3=1 3=1 3=1 

where it is required that a(e l9 ) = a(e~ ie ), and the tilde symbol denotes that only eigenvalues 
< 9j < it are considered in the average. 
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It remains to consider the cases of the orthogonal and unitary symplectic classical 
groups in (I2.10p . Their double scaling limit was first evaluated by Crescimanno et al. |19j . 
To present the results, it is convenient to introduce 

E N (L) oc £ ^< ■ ■ -^e-^^ 

ni,...,njv=— oo 

oo N 

E A 2 K...,^)(n-?)e- (w2/2L2)E - in ? (2.17) 

ni,...,n N = l j=l 

where the proportionality constants depend on N only Comparing with (I2.8P in the case 
A = 2n 2 /L 2 shows 

E N (L) OC e ^(N-l/2)(N-l)/^ L -2N^N z (^g^. +( 2 jV)) 

F N (L) OC e ^(N + l/2)(N + l)/6L^ L -2N^N z (^g^. Sp ( 2 jV)). (2.18) 

In terms of these quantities, the working of [19], further refined by calculations in [30] and 
[63] . gives that 

Urn ^log^(v^V(l +t/(2 7 / 3 iV 2 / 3 ))) = - l -{ql{t) -q' (t)), 

l^^hgF^Vmil+t/i^N^ = - l -{ql{t) + q' (t)). (2.19) 

Thus we have that (I2.10p the double scaling limits (I2.10p and (12.131) are identical, as are 
the results (I2.14p and (I2.15P to (I2.19p . Furthermore, as for the group integrals, only two 
of the three double scaling limits are independent due to the identity [26] 

oo 2N oo N 

e rufa) n k-^) 2 «( e rufa) n k 2 -^) 2 ) 

ni,...,n 2 ]v=— oo 1=1 l<j<k<2N ni,...,n N =—oo 1=1 l<j<k<N 

oo N 

x( e n k 2 -- 2 ) 2 )' ( 2 - 2 o) 

ni,...,njv=- oo 1=1 ^<j<k<N 

valid for general g even. 

The limit formulas reviewed in this section were derived as formally exact results, but 
without rigor. In distinction, rigorous working accompanied the limit formulas obtained 
in the setting of the next two subsections. 

2.3 Fluctuation formulas 

We now turn our attention to a seemingly unrelated class of problems, coming from sta- 
tistical physics. To introduce these problems, consider a unit square containing n points 
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with Poisson distribution Pr(n) = ^-e~ x * . Starting at (0,0), join dots with line segments 
of positive slope such that a continuous path ending at (1,1) is formed, and record the 
number of dots in the path. Call the maximum number of dots in such a path h u . This 
is a random variable which defines the so called Hammersley process (see e.g. J27J §10-9]). 
From work of Gessel [39] and Rains [60] we have that 



Pr(h D < N) = e- x2 (e XTr V +uf A (2.21) 

(the first of these references actually gives a Toeplitz integral formula, which can be shown 
to be equivalent to (I2.2ip ). This then provides a probabilistic interpretation of the single 
plaquette partition function in ( 12 .5p with gauge group equal to U(A r ). A rigorous analysis 
of Baik et al. [I] established that 

/ h n — 2X \ / f°° \ 

HmPrf ai/3 <tj =exp( - J (s - t)q 2 (s) dsj , (2.22) 

which one sees is completely consistent with (12.101) . 

Rains [60] formulated several symmetrized versions of the Hammersely process, in par- 
ticular by requiring that the points be symmetrical about the diagonal (line joining (0,0) 
to (1,1)), or the antidiagonal (line joining (0,1) to (1,0)). Again one wants to record the 
maximum number of dots in a directed path from (0,0) to (1,1). Let this number be 
denoted h m and h a respectively. Then we have from [60J that 

Pr(/i s <2A0 = e" A2 (e ATr5 \ 



/ SeSp(27V) 

Pr(fc <N) = V A2 ( (e ATr °\ + (e ATr °\ ) , (2.23) 

v ; 2 V\ /060+(7V) \ /oeo-(jv)/' v ; 

so providing a probabilistic interpretation of the single plaquette partition function in (12. 5 j) 
with gauge group equal to Sp(2iV) or O(N). Moreover, rigorous working in [6] established 
that 

aSL Pr (~W^ < ^ = 6XP ( ~ \ l t ( x ~ t ^ q2 ^ dx + \ f t <Io(z)dx^, 
, „ fh a -2X \ If / 1 f 00 , . 2/ , , 1 f°° . . , 
jimPrl A1/3 <t ) = 2\ eXp V 2 J ^ ~ ^ % ^ + 2 J 9o ^ 



+ 6XP (~^/ ( X ~ dx ~ \ J Qo(x)dx,^, (2.24) 

which, like the relation between ( 12.22p and (I2.10p . is completely consistent with (I2.19p . 

It is furthermore the case that the Yang-Mills partition functions ( 12. 8ft have interpre- 
tations in statistical physics [30], [63]. The most significant arises in the case of the gauge 
group Sp(2iV). A Brownian motion confined to the half-line x > and conditioned to start 
at x = when t = and to finish at x — 1 when t = 1 is referred to as a Brownian excur- 
sion. Suppose N such Brownian motions are placed at x = when t — 0, and furthermore 
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are conditioned to not intersect for < t < 1, before all returning to x = when t — 1. 
Let XAr(t) denote the position of the rightmost Brownian motion in [0, 1], and define the 
random variable = max{xN(t), < t < 1} specifying the maximum displacement. It 
was shown in [Ml Wf\ [22] that 

Pr (H N <L) = F N (L), (2.25) 
where F N (L) is specified by (I2.17P with the particular proportionality constant 

n 2N 2 +N 

An = ^ N2+N ' 2 UfJo r (2 + j)r(3/2 + j) ' 

Using the orthogonal polynomial method of [ID] and [19], it was shown by Forrester et 
al. [30] that 

lim F 7V (v / 2iV + 2^ 11/6 iV- 1/6 s) = exp ( -- [ (x - s)q 2 (x) dx + - [ q (x)dx). (2.26) 

at^oo \ 2j s 2 j s y 

Subsequently Liechty [H] used a Riemann-Hilbert approach to the orthogonal polynomial 
method to make this rigorous (see too the related work [T8]). 

The growth models reviewed in this section all have the feature that they give rise to 
expressions that appeared earlier in two-dimensional QCD. On should remark that there 
is now a vast literature on exact fluctuation formulas in growth models, particularly due to 
their relevance to the Kardar-Parisi-Zhang (KPZ) universality class (see [17] for a recent 
review) . 



2.4 Soft edge random matrix distributions 

The distribution functions exhibited in f l2.23f) and (I2.24p occur in a probabilistic setting 
again seemingly unrelated to what has been discussed to date, namely as the distribution 
of the extreme eigenvalue of certain universality classes of random matrices in the neigh- 
bourhood of a soft edge. By way of background, we remark that an Hermitian random 
matrix ensemble is said to have orthogonal (/3 = 1), unitary (/3 = 2) and unitary symplectic 
((3 = 4) symmetry if it has an eigenvalue probability density function of the form 

1 N 

T^n^) n i a *- a /- ( 2 - 27 ) 

N (=1 l<j<k<N 

The soft edge refers to a spectrum edge such that the leading order eigenvalue density 
decays to zero like a square root. As a concrete example, consider the case g(x) = e _/3x I 2 . 
Then (12.271) can be realized by Hermitian matrices with real (/3 = 1), complex (/3 = 2) 
and real quaternion ((3 = 4) entries (see e.g. [27J Ch. 1]). The entries on and above 
the diagonal are independent Guassians of mean zero and appropriate variance (which is 
different for the diagonal and off diagonal entries) . To leading order the eigenvalue density 
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is supported on [—y/2N, \/2N\, and vanishes like a square root at either end. Moreover, in 
the neighbourhood of such soft edges a well defined statistical state results by introducing 
the scaled variable s according to [2S] A = y/2N + s/{2 1 / 2 N 1 / & ). In particular the largest 
eigenvalue has a well defined limiting distribution in terms of these variables. Explicitly, 
let E^a^O; (a, b); g(X)) denote the probability that for the ensemble specified by (12.271) that 
are no eigenvalues in the interval (a, b). Then the distribution of the largest eigenvalue at 
the soft edge is specified by the limit 

£l oft (0; (s, oo)) := lim E PtN (0; (V2N + s/(2 1 / 2 N 1 / 6 ), oo); e"^ 2 / 2 ). (2.28) 

Beginning with a Fredholm determinant formula for (12. 28ft in the case /3 = 2 [25], Tracy 
and Widom [65] deduced that 



£| oft (0; (s, oo)) = exp ( - £ (x - s)q%{x) dx 



(2.29) 



(see [2], [MJ [12] for subsequent derivations). This result was known before (I2.22p . thus 
allowing for the interpretation of the latter in terms of the scaled largest eigenvalue of a 
complex Hermitian random matrix. Subsequently Tracy and Widom [66] gave Painleve 
formulas for the evaluation of the limit (I2.28P in the cases (3 = 1 and 4. The former of 
these reads 

E 1 soft (0;(s,oo)) = exp^-^ (x - s)q 2 (x) dx + q {x)dxj. (2.30) 



It was pointed out by Forrester and Rains [32] that the case (5 = A can be deduced from 
knowledge of the /3 = 1 and 2 results, and furthermore is most naturally expressed in terms 
of the rescaled quantity 

Ef t (0;(s,oc))= lim E p , N/2 (0; (V2N + S /(2 1 / 2 iV 1 / 6 ), oo); e"^ 2 ), (2.31) 
when it reads 

1 ( / 1 f°° If 00 \ 

Ef ft (0; (s, oo)) = - j exp ^ - - J (x - s)q 2 (x) dx + - J q (x) dxj 

+ ex p{~\J s ( x ~ s )^o( x ) dx - \ qo(x)dx^\. (2.32) 

Again, the evaluations (I2.30P and (I2.32|) allowed the results (I2.24p to be interpreted in 
terms of the scaled largest eigenvalue of a real, and real quaternion respectively, Hermitian 
random matrix. 

We remark that the matrix integrals in (12. 5p for each of the classical groups satisfy 
transformation identities relating them to the gap probabilities E^ ard (0; (0, s 2 ); a) with 
a = N and s proportional to bN [13]. The general /3 hard edge state with parameter a 
is realized by choosing g(X) = X a e~ l3X ' 2 , A > in (I2.27p . then scaling all the eigenvalues 
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A; H> Xi/(4N) and taking the limit N — > oo. It is called a hard edge state because the 
eigenvalue density is strictly zero for A < 0, and diverges like A -1 / 2 as A — > + . As a — > oo 
there is a hard-to-soft edge transition which gives an explanation for (I2.22p and (|2.24j) . 

A Painleve evaluation of E% (0; (s, oo)) seemingly different to ( 12 . 29[) was given by 
Forrester and Witte [31]. For this, introduce the Jimbo-Miwa-Okamoto cx-form of the 
Painleve II equation 

«) 2 + <(«T - K + u a ) - a 2 = 0. (2.33) 

Then it was shown that 

£f ft (0; (s,oo)) = exp ( - ^ u„(x) <fc) 

= exp ( — y ( s — x ) u o( x ) 5 (2.34) 

where uo(x) satisfies f !2.33|) with a = subject to the boundary condition 

u (x) ~ (Ai'(x)) 2 -x(Ai(x)) 2 . (2.35) 

X— ¥oo 

For future reference we remark that f!2.35|) implies that u' (x) has the asymptotic behaviour 

u' (x) ~ -(Ai(x)) 2 . (2.36) 

Theory relating to fl2T33|) [54] tells us that 

u' (t) = -2- 1/3 U /2 (t) + ql /2 (t) + |) (2.37) 



(see in particular [MJ eq. (5.27)]; note that in this work (II. ip is written with a h-> —a) where 
q a denotes a solution of the Painleve equation ( II. ip . and so we have that i?| oft (0; (s, oo)) 
can be evaluated in terms of go as specified by ( I2.29p . or in terms of q\/2 as specified by 
( I2.34p and ( I2.37p . As noted in [31], these two forms can be reconciled by invoking an 
identity of Gambier [38] which gives 

e2 1 /3 g 2(_ 2 -V3 t) = | g£/2 ( t) + eq 2 /2(t) + £ t> (2 . 38) 
valid for both e = ±1, and thus in particular u'(t; 0) = — <Zo(t). 



2.5 Parameter dependent problems and Lax pairs 

Underlying the analysis relating the Hammersely process to random matrix averages is 
a bijection between permutations and standard tableaux due to Robinson and Schensted 
(see e.g. [37]). As also detailed in [37], this was generalized by Knuth to give a bijection 
between weighted non-negative integer matrices and weighted semi-standard tableaux (or 
equivalently non-intersecting lattice paths). The latter underlies a class of lattice extensions 
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of the Hammersely model introduced by Johansson [45] and further generalized by Baik 
and Rains [TJ. These are defined by an ri\ x n 2 rectangular grid of lattice points 
each specifying a non- negative random variable Xij. The primary quantity of interest is 
the distribution of the maximum of the random variables when summed over an up/right 
path starting at the bottom left and finishing at the top right, referred to as the last 
passage time (this has the alternative interpretation as the maximum displacement in the 
non-intersecting lattice paths picture). 

In the case of a n x (n + 1) grid, Xij (j ^ n + 1) Poisson distributed with unit mean, 
Xi >n+ i Poisson distributed with mean (1 + ai)/2, the last passage time is distributed as for 
x\ in the probability density function proportional to 



- £?=i(*i+lH)/2-ai E"=i(zj-%)/2 



l<i<j<n 



(2.39) 



where x\ > yi > ■ ■ ■ > x n > y n > 0. It was proved in [5] that 



lim Pr (xi < 4n + 2(2n) 1/3 s 



f(s;w)Ef\0; (s,oo)), 



ai=-w/(2(2n) 1 / 3 ) 

where / relates to the Painleve II transcendent qo(s) via the Lax pair equations 



(2.40) 



d 


7(* 


w) 




ds 




w) 




d 


'f(s 


w) 




dw 




w) 





q (s) 
q (s) -w 



f(s;w) 
g(s;w) 



-wq [s J 



f(s;w) 
g(s;w) 



(2.41) 



wq {s) + q (s) w z 

(By definition of a Lax pair, qo(s) as it appear is (12.411) is determined by consistency of 

o^fo, and this gives rise to ( 11. ip with a = 0.) 



the mixed derivatives -t^t- and 



dsdw 



It is also the case that both / and g from the Lax pair appear as the distribution 
function for a particular symmetrized version of the above last passage time model [HI [7J. 
This is specified by a 2n x 2n grid with Xij (i < j) Poisson distributed with unit mean, 
for i > j and x^i Poisson distributed with mean (1 + a\)/2. The last passage 
time is distributed as for x\ in the probability density function proportional to 



l<«<J<2n 



(2.42) 



where X\ > x 2 > ■ ■ ■ > x 2n > 0, and it was proved in [6] that the distribution function 
permits the scaling limit 



lim Pr (xi <4n + 2(2n) 1/3 s 



(f + g)e- 



ai=-a/(2(2n) 1 / 3 ) 



+ (/ - g)e~* 



(Ef'iO; (s; (s, ex))))) 17 '. (2.43) 
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A random matrix interpretation of (I2.39P and (I2.42p . and thus of (I2.40P and (12.431) . 
was provided in [33]. Specifically, we have that (I2.39P with ai = 2/b — 1 is the eigenvalue 
probability density function of the rank 1 perturbed complex Wishart matrix, 

M = X ] X + bxx\ (2.44) 

where X (x) is an n x n (n x 1) matrix with standard complex Gaussian entries. With 
regards to (I2.42p . let X in (I2.44p be a 2n x 2n antisymmetric standard complex Gaussian 
matrix, and let x be a 2n x 1 standard complex Gaussian vector. Then we have that the 
perturbed eigenvalues of M have PDF (I2.42p with ai = 2/b — 1. 

It has been known since the work of Jones et al. [3S] that rank 1 perturbed Gaussian 
ensembles can exhibit a separation of the largest eigenvalue when the strength of the 
coupling exceeds a critical value. The significance of the coupling implied by (I2.40p and 
(I2.43P is that it corresponds to a scaling about this critical value, while the scaling of x\ 
centres about the leading order value of the largest eigenvalue, and measures in units such 
that the largest eigenvalues have spacing of order unity. The description of the rank 1 
perturbed scaled soft edge state for /3-ensembles ( (12.401) corresponds to (3 = 2; (12.431) to 
— 4) has been the subject of a number of recent papers [TO], [671 [281 El]- 

Another setting which gives rise to a Lax pair for qo{s) is the analogue of the double 
scaling limit for the so-called Hermitian quartic matrix model ( (12 . 2 7j) in the case /3 = 2, 
g(X) = e -(*/2)A 2 +( 9 /4)A 4 with ^ _ £2/(4^^2/3 = 0(1)) [g] . So does the scaled distribution 

of the maximum height in the nonintersecting Brownian excursion model detailed below 
(I2.24[) when considered also as a function of the scaled location of this maximum about 

t = 1/2 [621 ES GS]. 

We conclude this section by summarising a very recent calculation [69] of the joint 
distribution p s ^^(xi, x 2 ), of the scaled largest (xi) and second largest (x 2 ) eigenvalue at 
the soft edge for (3 = 2. To present this result, let p^f^it) denote the distribution of the 
scaled largest eigenvalue, so that 

P s 2 ?S)W = |^ oft (0;(t,oo)). 
Also, with q a a solution of (ll.ip . let 

Pa = q'a + Ql + \^ (2-45) 
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and introduce U(x, t) and V(x, t) through the Lax pair equations 



d_ 

dx 




-i 



x + 



-?3/2 



P3/2 



-1 



1 2 2 

- P3/2) + (?3/2 + ) 2 <?3/2 + — - 

2 £>3/2 P3/2 



1 P3/2 




'u 


-1 




V 



d 


~U~ 




"0 


0" 


di 


V 


- 


1 

.2 






X + 





-2(g 3 /j 



1 



P3/2' 





"[/" 


) 





(2.46) 



(it is noted in [69j Remark 1] that ( I2.46P are a transformed version of a Lax pair for 
PII first given by Flashka and Newell [23]; furthermore, as in [M], (11.11) is written with 
a 1 — y —a). With this notation, we have from j69j Prop. 9 and 10] that, subject to some 
specific boundary conditions for the transcendents involved, 



P2?2)(M - x) ~p s 2 %(t)t-^e-l t3/ \Ud x V - Vd x U){-2^ X] -2 1 / 3 *) 



An 
x exp 



dy (2g 3 /2 H 

2V3* \ P3/2 



)(-y) 



2y 



2?y 



(2.47) 



In the Appendix of [69] it is shown how q-s/2, through a combination of the (inverse) 
Gambier identity, 



g e/2 (-2 1 / 3 S ;0 = - e2 - 1 /3^ loggo(s;0; e = ±lj 



(2.48) 



and a Schlesinger transformation, can be related to the Hasting-Macleod solution q (see 
also Section [3T21 below) . 



3 Painleve II with general a in random matrix theory 
3.1 Two-dimensional QCD with quarks and matrix models 



In 1991 Minnahan [50] pointed out that the generalization of the one plaquette partition 
function (I2.5P with G = U(iV) for two-dimensional QCD without quarks to the case that 
each plaquette contains M massless quark flavours is given by 



» M 

I dUdtpdx exp (feiVTr (U + rf) + ^ x](I - XJ% + $(I - U ] )xi 
J i=i 



(3.1) 



where ipi and Xi an d N x 1 vectors with components which either commute (bosonic 
case, e = 1), or anticommute (fermionic case, e = —1). Moreover, it is noted that after 



13 



integrating out over ijj and Xi (13.11) reduces to 



exp (bNTi {U + U ] )t 2eM log | det(I - U)\) ) (3.2) 

V / / U&U(N) 

(of course as written (13 .2\\ is singular in the bosonic case and requires regularization). 

The quantity of interest is again the limiting "specific heat" as on the LHS of (12.101) . 
The method of orthogonal polynomials is applied in [50] to deduce that this evaluates 
in terms of a general Painleve II transcendent to give — It was stated that a is 

proportional to —M, although the proportionality was not determined. 

In 1998 Akemann et al. [3] studied the correlation kernel for the random matrix ensemble 
specified by (12.271) is the case g(\) = \\\ a e~ NV( - x \ where V(\) is such that the leading 
order eigenvalue density vanishes like a quadratic at A = (see [56l §13.2] for a text book 
treatment of the corresponding specific form of V). The correlation kernel (see (13.261) 
below, or more explicitly the orthonormal function entering therein) was shown to be 
given as the solution of a differential equation in which the coefficients are solutions of the 
general Painleve II transcendent with a proportional to a (the form of the general Painleve 
II equation in [31 eq. (2.50)] is not in the canonical form (II. ip . but doing so would give 
the explicit proportionality). Later Claeys et al. [15] gave a rigorous treatment of this 
problem using Riemann-Hilbert methods, and expressed the limiting correlation kernel in 
terms of the two-components of the Lax pair formulation of general a Painleve II due to 
Flascha and Newell [23] (this also came after the rigorous work of Bleher and Its [9] using 
Riemann-Hilbert methods to treat the case a = 0). Moreover, the Painleve II transcendent 
q a (s) appearing therein was specified to have the asymptotics 

q a (s) ~ J— - and q a (s) ~ - (3.3) 

s— >— oo y Z s— >-oo s 

and this solution was proved to be free of poles on the real axis, thus generalizing properties 
of the Hasting-Macleod a = solution appearing in Section 2. 

3.2 Random matrix averages 

The eigenvalue PDF for the Gaussian unitary ensemble (GUE) corresponds to the case 
g(\) = e~ x2 of ( l2~2Tj) . Let £* := 1 - e^(l - 0- The study of the averages 

^(A; n := ( n£i(l - £X(-oo,A))(A - W 

n£i(l-£*X(-oo,A))|A-A^) (3.4) 



where x^P = 1 f° r £ J i Xj = otherwise was undertaken by Forrester and Witte [3U [36] 
(for a later treatment, see [H]). Moreover, the corresponding soft edge scaled quantity 



^logi^ oft (s;/i;0 := ^log (e^ A2 / 2 ^(A;/i;0 



A^v / 27V+s/v / 2A rl/6 
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was computed in [35] (see also [27J Ch. 8]) with the result 

d 



, s \ogE s 2 ott (s;^0 = M^Oi (3-5) 



where satisfies (I2.33f) with a — ji. 
From the definitions 

Ef t (s-^ = 0^ = l)=Ef\0; (s,oo)), 

where i£| oft (0; (s, oo)) is given by (12.281) . and that uo(s;£ = 1) = u(s). For \i — and 
general £, by expressing _E| oft (s;/i = 0;£) as a Fredholm determinant [25] (see Section [3731 
below) it follows that the boundary condition (I2.35P is to be extended by multiplying by 

c, 

u (x; ~ £((Ai'(z)) 2 - x(Ai(*)) 2 ) ~ -^e"^ 2 . (3.6) 
Furthermore, as a consequence of the Gambier identity (I2.38p . and (I2.36p . we have 

u' (t;0 = -q 2 o(t;0, (3-7) 
where qo(t;£) satisfies fll.ip with a = subject to the boundary condition 

-(2/3)i 3 / 2 

qo(t;C) ~ v / ?Ai(t) ~ v /£^^. (3.8) 

We remark that the quantity i?J (s; /x = 0; £), < £ < 1, has the probabilistic interpreta- 
tion that in the /3 = 2 soft edge state, with each of the eigenvalues deleted with probability 
(1 — £), there are no eigenvalues in (s, 00) [TTj . 

Note that (13.81) with £ = 1 is consistent with (12.111) . In this case we know that [32] 

»>(*;£ = 1) ~ (3.9) 
t— >— 00 y z 

and this used in (I2.29P has the consequence that i?2 Oft (0; (s, 00)) ~ e -(- s ) 3 / 12 |65j . On 

S— OO 

the other hand, it is known [1] that for < £ < 1, qo(t; £) decays to zero as £ — >■ — 00 and 
oscillates, 

gb(f;0, ~ dHr 1/4 sm(^H) 3/2 -^ 2 log(-t)-c) (3.10) 

where 

1 3 1 7T 

d 2 = — log(l - 0, c= -rf 2 log2 + ar g r(l - -id 2 ) - - 

7T Z Z 4 

and this implies the very different asymptotic behaviour |ll) . [2"9~| 

j* = 0;0 ~ exp(log(l-0(2/3vr)(- S ) 3 / 2 ). 
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We seek the boundary conditions relating to ( 13.51) for general £ and values of /i different 
to zero. In the cases \x = 1 and \x = 2 this can be done by using known formulas which 
give relations with the case \x = 0. In particular we require the (inverse) Gambier identity 
(I2.48p . together with the recurrence [3U Eqns. (3.23) and (5.9)] 



_x(x; - u^x- = -2 1 / 3 g At+1/2 (-2 1 / 3 a;; f), (3.11) 
and the identity [Ml Eq. (5.22)] 

«2(s;0 = ^: lo s n o(^;0 + «o(s;0- ( 3 - 12 ) 

Proposition 1. We have 

15 1 — (4/3)2.' 3//2 i 

and 

u 2 ( X ; Q ^ -2^ -- + — + (-^ + f (l + O(^) ) . (3.14) 

Proof. We set /i = in (13.111) . then substitute for q 1/2 (-2 1/3 x; £) using (12. 48ft with e = 1 
to deduce 

tti(x;0 = w (x;£) + — logg (s;0- ( 3 - 15 ) 
as 

Next, we extend (13. 6p by writing 

/ >\ £ 4 t3/2/, c Cl 1 



8vrx V 2x 3 / 2 2x 3 Vx 9 / 2 

+ tS^^ ( i+ i4 +o (i)) + °(^n < 3i6 > 

Substituting for a in (12.331) with a = shows 

17 /35\2 13 , 1531 , 01 _ 

c = , ci = — , a = , d\ = . (3.17) 

12 V24/ 6 288 ( J 

Furthermore, we know from [58] in the case £ = 1, and [271 Prop. 9.4.4]) for general £ > 
that the analogous extension of ( 13. 8 p is 

*( s;f ^~^v^i 1 -# +0 t) + ^6^( 1 -# +0 b))}- (318) 

where 

5 23 
ai =72' 01 = 24' 
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Substituting ( 13 . 1 8[) and ( 13 . 1 6[) in (I3.15P gives, after straightforward manipulation, the ex- 
pansion (13.131) . 

To obtain the expansion (13.141) . we only require (I3.16p . Thus we substitute in ( I3.12p 
and expand the logarithm to obtain the stated result. □ 

The structure exhibited in ( 13 . 1 3[) and ( I3.14p . where the dependence on the parameter £, 
which itself is not a parameter in the Painleve II equation (II. ip . occurs only in exponentially 
terms is well known in the theory of Painleve equations (see e.g. [21])) • The full asymptotic 
series of such terms is said to form a trans-series. 

The leading x — > oo asymptotic expansion of u^(x;l;) for fi = ±1 and < £ < 1 can 
also be determined. In addition to (13. lip , this requires the expression for u(x; y«;£) (which 
is essentially the Hamiltonian for the PII system in Okamoto's theory [51]) in terms of 
q»(x;Z) [34j Eq. (5.26)] 

- 2- 1 /3 ?v (_2^/3 t . = ^( q ' 1/2 ^(t;0) 2 -\{qI/2-^0 + 1) 2 - (3.19) 

Proposition 2. We have 

x l/2 

u ±1/2 (x;£ = l) ~ T^-, (3.20) 

while for < ( < 1 we have 

-2-V3 (2-V3 x . „ Ml ~ x l/2 CQS 2 ft 3/2 _ 3 d2 x _\ 

x-^oc 271 \3 4 / 

where d and c are as in A3.10\) . 

Proof. The behaviour of uyzfa; 1) given in ( I3.20p follows by substituting \i = 1/2 in 
( I3.19p . then using ( 13. 9p . With this established, the behaviour of tt_ 1 / 2 (x; 1) follows after 
noting from (13. lip that 

u 1/2 (x; - u_ 1/2 (x; = -2 1 / 3 g (-2 1 /3 x; ^) ( 3 .22) 

and a further use of (13. 9p . 

To derive (13.211) in the case fi = 1/2, we simply make use of (13. 19[) with fi = 1/2, and 
then substitute (13. 10p . The result for /i = —1/2 now follows from (13.221) and further use 
of (13T0D . □ 

In distinction to the forms ( 13 . 1 3f) and ( 13. 14}) for fi = 1 and 2, the dependence on £ 
for /i = ±1/2 as exhibited in Proposition [2] appears already in the leading term, and 
furthermore this term is, for £ ^ 1, oscillatory. This latter feature can be understood from 
the details of the average (13. 4p . Thus for £ = 1 the quantity being averaged is always 
positive independent of /z, while for < £ < 1 and \x = ±1/2 this quantity can take on 
complex values, but for /i G Z + and £ in this range it is again real. This suggests that in 
the variable £* as defined above (13.41) . with < £* < 1, the leading asymptotic form will 
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be analogous to that seen in ( 13 . 1 3[) and (I3.14p . and that the £ dependence will only be seen 
in a trans-series. As already noted in [3U Eq. (5.35)], this leading form is expected to be 



u^{x]i) fix 1 - — + . +■■■ 3.23 

According to (13.191) and (13.221) this could be checked for fi = ±1/2 upon knowledge of the 
x — > — oo form of qo(x, 1 ± i(l — £*)), if this was to be available (see [TB] for an easy to 
read review on what is available). 

For general (but fixed) values of £, Prahofer and Spohn [58] have shown how the Hasting- 
MacLeod solution qo(s; £) can be numerically computed to very high precision. According 
to formulas presented above, this means too that u M (x;^) for fi = ±1/2,1 and 2 can 
similarly be numerically computed. Moreover, with these values as initial conditions, we 
can then compute u^(x;^) for general positive integer or half integer values by using the 
fact that the latter satisfy the alternate discrete Painleve I equation [3H Eq. (5.12)] 

+ 1 1 < . * = * " ( W*; - 0) 2 . (3.24) 



u^+i(x; f ) - u^i(x; Cl u^X] f) - u^x; 

We remark too that this equation, in theory at least, allows the leading ^-dependent terms 
in the x — > oo asymptotic expansion of n M (x; ^) to be determined for all \\i = 0, 1, 2, ... . 
However in practice, as it appears from (13. 6p . (13.131) and (13. 14|) that the leading such term 
is proportional to ££" 1-3M / 2 e~3 x3/2 ; this effectively means expanding the trans-series of 
Uq(x; £) to one higher order for each integer increase in /i, due to cancellations. 



3.3 Boundary spectrum singularity 

Closely related to (E3D is the PDF 11237]) in the case (3 = 2 and g(X) = |A - a^e'^ . 
For fi G Z + , this has the interpretation as the eigenvalue PDF for the Gaussian unitary 
ensemble, conditioned so that the point A = a is an eigenvalue with multiplicity /i. It is 
well defined for general fi > — 1. 

The best known case of this PDF is when a = 0. In [53] the correlation functions for the 
scaled state defined by replacing each A/ by nXi/V^N and taking iV — > oo were computed. 
In keeping with the interpretation for \i £ Z + , this state was subsequently referred to as 
the (bulk) spectrum singularity (see e.g. [27]). More recently, Its et. al [13] considered 
the situation in which the spectrum singularity at A = a is located within the soft edge 
boundary layer. (In fact in [13] the generalization of PDF with e~ x replaced by e~ v<yX \ 
for any V(x) inducing a soft edge, was considered.) 

In general the /c-point correlation function for a unitary invariant ensemble (/3 = 2 case 
of (12.271) ) is given by a determinant 

p( k )(xi, ■ ■ ■ ,x k ) = <tet[K N (xj,xi)] jtl=lt ... t k, (3.25) 
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where the correlation kernel is given in terms of the orthogonal polynomials of one variable 
{Pn( x )}n=o,i,- w ith respect to the weight function g(x) according to 

N-l 1 

K N (x,y) = {g{x)g(y)) l/2 ^ TFPj( x )Pj(v) 

= (9(x)g(y))^ PN{x)pN ^ y) 7 p»Mpn-i{x) 

Mn-i(x - y) 

(see e.g. [27J Ch. 5]) with Afj := g(x)(pj(x)) 2 dx. The main result of [13] is the 
evaluation of the scaled limit 



K s u oit (x,y;c) := lim * - K N (V2N + J* - , v / 2iV+^- , 
The functional form 



/2iV+c/(\/2V 1 /6) 



tp 2 (x\ c)^i(y; c) - ViO; c)^ 2 (y; c) 
(X ' y;Cj " 27ri(x-y) [6 - Z7) 



was obtained, where the ipj, given explicitly in [13] . can be expressed in terms of the 
Flaschka-Newell [23] Lax pair components for the Painleve II equation (II. ip with a = 
ft + 1/2. 

This characterization was shown to dramatically simplify in the cases \x = 0, 2 to forms 
involving only Airy functions, 

^soft, „,. x _ Ai(x + c)Ai'(? / + c) - Ai'(x + c)Ai(j/ + c) 



x-y 

^soft 



K% oft (x, y; c) J^ oft (0, 0; c) - i^ oft (x, 0; c)i^ oft (?/, 0; c) 
^ oft (0,0;c) 



These simplification were achieved by firstly demonstrating that the Flaschka-Newell Lax 
pair can be written in terms of Airy functions for a = 1/2. 

The functional form for KQ 0it (x,y; c) is the fundamental Airy kernel (in the variables 
x + c, y + c)) which underlies (I2.29p . The corresponding fc-point correlation function 
p\*k^° * s gi ven by substituting this for Kn in (13.251) . On the other hand, we see from the 
definitions that for m — 0, 1, . . . 

(2m+2,c)/ \ 
(2m+2,c), \ i • ^(fc+1) 1^1) ■ ■ -^fej^fe+l) ; /o on\ 

P fc) ; (x 1 ,...,x fc )= hm ^— (3.29) 

Xfc+1 ^ U P( 1} '(Xfc+l) 

(the case m = of this formula was used in [31] in relation to the bulk spectrum singularity 
ensemble). This allows us to express -Klm+2 m terms of thus generalizing the second 

of the formulas in (13. 28j) . 
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Proposition 3. We have 



K^ +2 (x,y; c) = lim 



K^(x, y- c)K%?(z, z- c) - K?«(x, z; c)K^(z, y; c) 



z->0 



K^(z,z;c) 



(3.30) 



Proof. We have, for small Xk+i, 

(2m+2,c)/ \ 
P(k+1) \ X 1: ■ ■ ■ : X k,Xk+l,C) 



-- det 
det 



[K%£(Xj, Xu c)]j, 1=1,...,*. [i^2m(^, 0; c)] j=lj ... tk 

_ [K^(0,x k ;c)]i = i y .. tk K 2m {x k+1 ,x k+1 ;c) 

[ K ?m+2{Xi,Xh C)]j,l=l,..., k [0]j=l,...,k 

[0]«=i,..., fc K%£(x k+1 ,x k+1 ;c) 



where the second equality follows by applying elementary row and column operations to 
the first determinant., and taking x k+ i — » in the first block. Expanding by the last row, 
and noting that 

P(i)"' c) (^fc+i) = K2^(x k+1 ,x k+1 ;c) 

gives (1X50]) . □ 

It is a basic but fundamental fact that the gap probability in a statistical mechanical 
system can be expressed as a sum over the correlation functions (see e.g. [27J Eq. (9.4)]). For 
the boundary spectrum singularity ensemble, the generating function for the probability 
of exactly k (k = 0, 1, 2, . . . ) eigenvalues in the interval (0, oo) is precisely the ratio 



£S0ft ( 



£! oft (c;/i;£ = 0) 



where i^ oft (c; /x; £) is as specified in Section [3.31 Using the expansion of this generating 
function in terms of the correlations one has 



£S0ft ( 



s;//;£) 



£f ft (s;/i;£ = 0) 
Recalling (13. 5p . it follows that 



(3.3i; 



Wm(c;0 - w M (c;^ = 0) ~ -£ 

c— 5>oo 



^x) dx. 



(3.32) 



In the case \i = 0, (13.321) is equivalent to (13. 6p . Moreover, we can use f!3.32j) to re- 
claim the £ dependent term in (13.141) . Thus a straightforward calculation (peformed using 
computer algebra), expanding K^ iXl {x,y] c) as specified by (13.281) for large c shows that 



K, 



)ft (x,y;c) ~ e -W^ 3/2 e -(-+y)V-cjy_ 

c— >oo 2 TIC 
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Recalling that p%. (x) = i^| oft (x, y; c) we substitute this in (I3.32p to conclude 



e -(A/3)c^ 
c^oo 13 2567TC 4 ' 



M c ;0 - "m( c ;£ = 0) o CC _^ > ( 3 - 33 ) 



which is in precise agreement with (I3.14p . 

We can can furthermore demonstrate that all terms proportional to £ as given by 
f)3.32p agree with those obtained from (I3.12p . For this we must extend the first asymptotic 
equality in (13. 6p to read 



«o(c;0 



POO 

~ ^ oft (0, 0; c) - £ 2 / (^ oft (x, 0; c)) 2 rfx + 0(£ 3 ). (3-34) 

This is a consequence of (13.51) and the fact that i?| oft (s; 0; £) = det(l — £-K( Sj00 ), where 
-K"(s,oo) is the integral operator on (s, 00) with kernel i^Q° ft (x, y; 0), and the expansion of 
the latter in terms of determinants (or equivalently the underlying correlation function; 
essentially we are extending the case /x = of (I3.3ip to second order — as an aside we 
remark that this formula provides an alternative method to derive the coefficients in the 
expansion (13.161) ). Substituting (13.341) in (I3.12p reclaims, after minor manipulation, (13.32)) . 
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